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Low and high energy properties of the pion are reviewed in the frame- 
work of chiral quark models. Particular emphasis is put on the simplest 
' version of the SU(2) NJL model as prototype. The role of gauge invari- 

ance in this kind of calculations is stressed. The results are used as initial 
CN ■ conditions for perturbative QCD evolution equations. At leading order the 

' quark model scale is /zq 320MeV as determined from the pion distribution 

functions and the pion distribution amplitudes. 
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1. Introduction and Motivation 



The structure and dynamics of hadrons is known to be a difficult problem 
because tiie underlying QCD theory in terms of quark and gluon degrees of 
\ freedom is actually very complicated [Q]. There are, however, well known 

situations where some helpful simplifications take place. 

At high energies asymptotic freedom sets in, and perturbative QCD 
becomes applicable. The renormalization group equations allow, on the 
basis of the operator product expansion (OPE), to relate hadronic matrix 
elements of quark and gluon operators at different scales fH . There are many 
situations where this scale dependence has been predicted and in some cases 
tested, particularly in the study of high energy inclusive Q and exclusive 
processes Q, but perturbative QCD says nothing about the value of these 
matrix elements at a given scale. 

At very low energies, spontaneous chiral symmetry breaking of the QCD 
Lagrangian dominates, and the would-be Goldstone bosons in the massless 
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quark limit become the lightest states in the hadronic spectrum. They are 
naturally identified with the pion, and it can be shown that at low energies 
the interaction among pions becomes small. Actually, at energies of the 
order of the pion mass, Chiral Perturbation Theory (ChPT) can be applied, 
in terms of low energy constants which correspond to certain QCD Green's 
functions, and many successful results have been found |5|, ^. The values 
of the low energy constants cannot be computed themselves from ChPT. 

The previous considerations suggest the study of pion properties at high 
energies, be exclusive or inclusive processes, trying to exploit Chiral Symme- 
try. From a theoretical viewpoint, vr mesons are particularly distinguished 
hadrons since most of their low energy properties follow the patterns dic- 
tated by chiral symmetry. Actually, we do not expect to understand the 
properties of any hadron better than the pion, as Chiral Perturbation The- 
ory suggests. Recent studies have shown, for instance, that ChPT can be 
applied to compute chiral corrections of hadronic matrix elements involved 
in deep inelastic scattering Q (in much the same way as chiral corrections 
are computed for vrvr scattering). However, nothing is said on the value of 
those matrix elements in the chiral limit. 

The calculation of pion matrix elements of quark and gluon operators 
involved in high energy exclusive and inclusive processes requires some non- 
perturbative method. Along these lectures it will be shown how this can 
be done with the help of low energy chiral quark models, which account 
for spontaneous chiral symmetry breaking, supplemented with QCD evolu- 
tion equations, which automatically implement QCD logarithmic radiative 
corrections. 

Chiral Quark Models share with QCD several features, mainly chiral 
symmetry (realized in the Goldstone phase) and explicit quark degrees of 
freedom. As we will see below they are useful in that they can make quanti- 
tative predictions which, in many cases, agree rather well with experiment. 
Unfortunately, there is not an unique model, but many and almost as many 
versions as authors working on the field. So, for practical reasons we will be 
using the Nambu-Jona-Lasinio model |^] for quark degrees of freedom 
It turns out that, when dealing with low energy properties of the pion, 
model details seem rather irrelevant, so most people have not paid partic- 
ular attention to their different versions. This seems not to be the case 
for properties which are extracted from high energy inclusive and exclusive 
processes where the model details and/or their different versions turn out 
to produce quite different results. 

Given the fact that QCD is a theory with quarks and gluons, one obvi- 
ous question one should ask is what is the scale at which a (chiral) quark 
model is naturally defined taking into account that its only explicit degrees 
of freedom are quarks. So far, we only know of one possible quantitative an- 
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swer [10, namely the adequate scale is defined by requiring that quarks 
carry all momentum of the hadron. In a field theoretical chiral quark model 
this is a trivial statement, because there are no explicit gluons. This as- 
sumption seems consistent with QCD in perturbation theory since it turns 
out that the gluon contribution to the energy-momentum tensor becomes 
smaller as the running scale goes down until it eventually vanishes. The 
remaining quark contribution can be separated into sea (flavor singlet) and 
valence (flavor non-singlet) pieces. The precise value suggested by pertur- 
bative calculations is in the range, /i = 300 — 350MeV. Of course, the 
separation of the gluon, sea-quark and valence-quark momentum fractions 
is a not a renormalization group invariant operation. Thus, if such a match- 
ing procedure between QCD and the chiral quark is meaningful, it depends 
on the particular renormalization scheme. 

2. Low energy models from high energies 

QCD is an asymptotically free theory, which means that the universal 
coupling constant among quarks and gluons becomes small as the renor- 
malization scale goes to infinity, and hence perturbation theory becomes 
applicable. For definiteness, we use LO evolution for the running strong 
coupling constant 

= (I) l0g(.v'A?,CD) 

with (3q = UCa/'^ — 2Np/3, Ca = 3 and Np being the number of active 
flavors, which we take equal to three. We take for concreteness Aqcd = 
226 MeV, which for n = 2GeV yields a = 0.32 Q. In perturbative QCD 



calculations, a/27r is the expansion parameter. Renormalization group in- 
variance requires a very specific dependence for any set of observables, O, 
on the scale 

0{fi) = Uiii,fio)Oipo) (2) 

where U{fi, /xq) is a linear matrix operator, fulfilling cocycle properties typ- 
ical of evolution equations, 

U{ni,fl2)U{fl2,fJ'3) = U{fli,H3) (3) 

C/(/x,/i) = 1 (4) 

In the limit /i' — > oo the operators are explicitly calculable in perturbation 
theory. Obviously, conserved quantities do not depend on the evolution 
scale. 
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A particularly interesting operator in QCD is the energy-momentum 
tensor 9°^^, which due to Poincare invariance is a conserved quantity and 
hence renormalization invariant. Its matrix element between pion states 
fulfins, 

(vrir^l^) = 2pV (5) 

For the QCD Lagrangian the energy momentum tensor can be separated 
into three contributions Q, 

B-^ = e-J + ef, + ef^, (6) 

where Oq^ ,9q^s and 6'°ns gluon, quark-singlet and quark- non-singlet 

(or valence) contributions respectively fulfilling 

{7r\efs\^)\, = 2pV5i(/i) (7) 

where Gi(/x), 5'i(^) and Vi(/i) are the gluon, sea quark and valence quark 
momentum fractions of the pion. In deep inelastic scattering (DIS) it can 
be shown Q that if qi(x,/i), qi{x,fi) and G{x,^) represent the probability 
density of finding a quark, antiquark and gluon respectively with momentum 
fraction x at the scale // (jj? = in DIS) then one has 

Giin) = f\xxG{x,fj,) (8) 
Jo 

'SiifJ') = / dxxS{x,iJ,) (9) 
Jo 

Viifi) = [ dxxV{x,n) (10) 
Jo 

where, for 7r+, V = Un — Un + dn — cItt and S = u-,^ — Un — dj^ + dj^. These 
momentum fractions depend on the scale, ^, and fulfill the momentum sum 
rule 

Viifi) + Siifi) + Giif,) = I (11) 

as a consequence of the energy-momentum tensor conservation. In per- 
turbation theory it is verified that due to radiative corrections Gi{fi) and 
Si{fi) decrease as the scale /i goes down. On the contrary, the non-singlet 
contribution to the energy momentum tensor evolves as 



^i(^o) \a(/^o) 
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where 7i /2/3o = 32/81 for Np = Nc = 3. The value of Vi{n) has been 
extracted from the analysis of high energy experiments. In Ref. [|l^ it was 
found that at = 4GeV^ valence quarks carry 47% of the total momentum 
fraction in the pion, e.g., for 7r+, 

Vi = {x{u^-u^ + d^-d^)) = 0A7 ±0.02 at /i^ = 4GeV2 . (13) 

Downwards LO evolution yields that for a given reference scale, f^o, 

Vi{fio) = l GM + Si{fio) = (14) 

The scale no so defined is called the quark model point for obvious reasons. 
At LO the scale turns out to be 

no = SlStio^eV (15) 

This is admittedly a rather low scale, but one can still hope that the typical 
expansion parameter a(//o)/2vr ~ 0.34 it 0.04 makes perturbation theory 
meaningful. Since such an approach was first suggested jl^ that is all 
one can do for the moment Q. There are more sources of uncertainties 
to Eq. (|T5|). For instance, if the point for the quark model is defined by 
Giinc) = 0, then at LO the scale is fic ~ 350MeV. The determination of 
Ho given by Eq. (^) is model independent. In these lectures we will show 
how this determination not only leads to a successful description of non- 
singlet pionic parton distribution functions in certain versions of the chiral 
quark model, but also that the number is in quantitative agreement with 
other determinations. 

3. Chiral symmetry and chiral quark models 

For the SU(2) up and down quarks, chiral symmetry is the invariance of 
the QCD Lagrangian under the transformations Q, 

q{x) e''^q{x) (16) 
q{x) ^ e'"-^g(x) (17) 
q{x) e'^-^^"q{x) (18) 

where q{x) represent Dirac spinor fields with Np = 2 flavors and Nc colors 
where r are the isospin Pauli matrices. As a consequence of this symmetry 

^ Actually, in the case of the nucleon these low scales produce negative gluon densities, 
if one takes Vi{n — 2GeV) — 0.40 and hence violate positivity of parton distribution 
functions as well as unitarity of structure functions |14l lla]. 
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there exist baryon, vector and axial Noether currents, 





q{x)'^^q{x) 


Baryon 


current 


(19) 






Vector 


current 


(20) 






Axial 


current 


(21) 



respectively. Conservation of the vector current (CVC) and partial conser- 
vation of the axial current (PCAC) imply that 

(22) 
(23) 

mq{x)i^5Taqix) (24) 

with m denoting the average current quark mass (we neglect isospin break- 
ing effects). A chiral quark model is any chirally invariant dynamical field 
theoretical model containing only explicit quark degrees of freedom and ful- 
filling the previous conservation laws ([2^), (p^) and (p^. Unfortunately, 
there is no such a thing as the chiral quark model. This is probably the rea- 
son why there exists a proliferation of models with many variants. Regard- 
less of their differences all these models are characterized by the three ways 
of breaking chiral symmetry in QCD: explicit, spontaneous and anomalous. 
These are non trivial requirements, which imply constraints on the regular- 
ization methods in local models or equivalently on the high energy behavior 
of non-local models. In fact non-local models have been preferred because 
they provide a more natural explanation of the anomalous vr'^ — > 77 de- 
cay rate ||l^, but the calculations are cumbersome. In addition, they are 
formulated in Euclidean space and their extrapolation to compute matrix 
elements involved in high energy processes is subtle as it is on the lattice. 
For instance, the calculation of structure functions requires either a contin- 
uation of the non-local model to Minkowki space or the determination of 
some moments in the Euclidean region and subsequent distribution recon- 
struction. While the first possibility turns out to be extremely difficult to 
explore, the second alternative cannot be used to pin down the x — > 1 which 
depends on the asymptotic behavior of the moments. On the contrary, local 
models, although less realistic, require specifying a suitable regularization 
which may be directly formulated in Minkowski space. 

3.1. NJL Lagrangian 

For the purpose of our discussion we think it is useful to review the 
SU(2) NJL model ^J. This model has extensively been used in the past 



df.Jy'^ix) = 
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and there exist many reviews on the subject. So, although we aim at a self- 
contained discussion of pion properties, our presentation will be necessarily 
sketchy in order to provide the main ideas. Nevertheless, we will stress those 
points where differences with other authors become important. The NJL 
Lagrangian in Minkowski space is given by ^ 

Cmi. = q{i9 -m)q+- ({qqf + {qfi^^qf) (25) 

where q{x) is a Dirac spinor field with N^. = 2> colors and Np = 2 fla- 
vors, G is the coupling constant with dimension [G\ = Af^^ and m the 
average mass of the current up and down quarks (~ 7MeV at the scale 
IX ~ IGeV). With the exception of the term in the current quark mass m 
, this Lagrangian is invariant under the UsiX) ® SU{2)r SU{2)l chiral 
group, with a U{1)b ® SU{2)y subgroup. Thus, the currents (pj|), ( pO[ ) 
and (^T|), and their conservation laws (p^), (|23|) and ( p^ are satisfied, as 
in QCD. A consequence of Poincare invariance is the conservation of the 
energy momentum tensor, 

e^'^' = ^q {j^d" + Yd") q - g^" C (26) 

which only contains quark degrees of freedom. Thus, we expect in this 
model that for any hadron, and in particular for the pion 

(vrir^lvr) = 2p>'^. (27) 

Since the coupling constant, G, has a negative energy dimension, the La- 



grangian (25) above is not renormalizable by power counting and it is usually 
interpreted as a theory with a finite cut-off, which we generically denote by 
A. This means that results depend crucially on the value of A and some fine 
tuning will be invoked. The need of a cut-off is not a serious problem from 
a physical point of view because we know that at high energies, much above 
the cut-off, the effective interaction should be replaced by the underlying 
QCD interactions in terms of explicit quarks and gluons. Roughly speak- 
ing, this feature is present in more sophisticated chiral quark models with 
non-local interactions. 

The real problem derived from the high-energy suppression has to do 
with the fact that we want to preserve the symmetries of the Lagrangian, 
namely gauge and chiral symmetry, without removing the cut-off by taking 
the limit A — > oo, i.e. keeping the high energy suppression. This implies 
that not every regularization, although it may accommodate some particular 
prejudices, may be considered as suitable. This point will be illustrated 
below. A full discussion on regularization methods in the NJL model can 



be found in Ref. 
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In these lectures we will work, as usual, in the one-quark loop approxima- 
tion. This is equivalent to work to leading order in the large Nc expansion, 
taking GNc =constant. In this way we comply not only with chiral symme- 
try constraints, but also with large Nc requirements. The study of higher 



orders in a 1/Nc expansion can be traced from Refs. |18, 19, [2 



3.2. Pauli-Villars regularization 

In a cut-off theory most of the discussion is related to the regularization, 
and hence gets a bit technical. We will use throughout these lectures the 



time honored Pauli-Villars (PV) regularization method |21] which has al- 



ready been used for several applications in the NJL model, like ChPT p3 



finite temperature ||2^ , solitons , structure functions |£5|, |2^ , correla- 



tion functions |26] and distribution amplitudes |2E]. To our knowledge this 
is the most reliable regularization method of the NJL model so far fulfill- 
ing many desirable features, among others preserving gauge invariance and 
working directly in Minkowski space (see the recent discussion in Ref. [p8| ). 
The only difference with the standard PV method of QED is that in the 
NJL model this type of regularization is applied at the constituent quark 
mass level (see below). It is surprising that given the advantages of the 
method, it has been used so few times as compared with other methods 
such as Euclidean 0(4) cut-off or Schwinger's Proper-Time. A potential 
disadvantage of the method is that it provides a negative spectral strength 
due to the PV subtractions, and thus positivity of some physical quantities, 
like form factors or distribution functions might not be fulfilled. Although 
this is not exclude, in all the calculations we present we explicitly see that 
is not the case. 

In the NJL model this method has been used mainly within the context 
of the bosonization or auxiliary field method (see e.g. Ref. ||2^), but at the 
one loop level and for the processes we will be considering it corresponds to 
making the replacement under the momentum integral 

1 1 , / h + M _Jn±M_\ 



]/tl-M If^N-M ^ '\kf-M^-A^ k%-M'^-Af 

Ci and Aj are the same suitable coefficients for all one loop graphs, fulfilling 
Co = 1 and Aq = and chosen in such a way as to make one loop integrals 
finite. In the graphs with an odd number of Dirac 75 's the regulator has to 
be removed at the end of the calculation. 

In the NJL there appear quadratic and logarithmic ultraviolet diver- 
gences. To illustrate the method let us consider the integral 

f d^k 1 
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which is quadraticahy ultraviolet divergent, as can be seen by integrating 
first over ko component (using Feynman's — > — iO"*" prescription) 
and then over the k component. Using the PV subtractions, it becomes 

^^^ = "^^^-/ (t)4fc^-M2-Af 

Evaluating the integral by any method, one sees that at least the following 
two conditions are required to render the one loop Feynman integrals finite 

Y.Ci = Eci^' = (31) 

i i 

Thus, at least two subtractions are needed. Solving the system of equations 
in terms of Ai and A2 provides two unknown parameters. To reduce them 
to only one the coincidence limit Ai ^ A2 = A is taken yielding the rule 

J2 c^fiAf) = /(O) - /(A2) + A2/'(A2) (32) 

i 

This calculation already illustrates a very general feature of the PV method, 
calculations may directly be done in Minkowski space, and there is no need 
to go to Euclidean space. This is a computational advantage in the study 
of high energy processes as we will see below. To simplify the notation we 
will implicitly assume the use of such a regulator in what follows. 



3.3. Chiral symmetry breaking 

To illustrate the method let us consider the Dyson-Schwinger equation 
for the quark propagator, S{p), (see Fig. |l]). 

S{p) = So{p) + S{p){-i) J ^Tr [S{k)G] So{p) (33) 

where S{p) and So{p) are the full and free quark propagators respec- 
tively. G is the four-point kernel, which in the lowest order approximation 
becomes the coupling constant and can be taken out of the integral. Thus 
the full propagator acquires a constant self-energy, S(p) = M. Computing 
the Dirac, flavor and color traces one gets the so-called gap equation in the 
convenient form 

M (1 - AGNfNJ) = m (34) 

where Eq. (pO|) is understood for the integral /. For m = there are two 
solutions (the Wigner alternative) 
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Fig. 1. Schwinger-Dyson equation for the quark propagator. The thick hne repre- 
sents the fuh propagator, whereas the thin Hne stands for the free quark propagator. 
The full blob is the irreducible two particle amplitude. 



• Wigner phase. Corresponding to M = 0. 

• Goldstone phase. M 7^ and hence 



1 = AGNpNJ (35) 

which is called the gap equation in the limit m —> ^ because it provides 
a gap of width 2M in the quark spectrum. 

For m 7^ the gap equation is Eq. ( |3^ ) and one has to choose the solution 
which goes to Eq. p5|) for m — > 0. 



3.4- Quark condensate 
In the chirally broken phase, one has a quark condensate ( {qq) = {uu + 



dd)) 



I , / x„ / d'^k 1 M — m, , ^ 



The standard accepted value in from QCD sum rules |Q is {qq)fj, = {uu + 



dd)f, = -2(240 ± lOMeV)'^ at the scale /i = IGeV. At leading order in the 
QCD evolution one has 

(g?)^ _ "(^) V''^^* 7- -8 (37) 

If we use the scale of Eq. (|15|) we get an enhancement factor ~ 2 and hence, 

= -(380±20MeV)l (38) 

According to our point of view it is this value the one that should be com- 
pared to the model calculation, and not to the one at /i = IGeV. Of course, 
one may argue that LO evolution is not sufficient to go to such low scales. 
All we can do to check the consistency of the approach is to provide, as we 
do below, other determinations of the scale fiQ. 
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Fig. 2. Bethe-Salpeter equation for quark-antiquark scattering. The square blob 
represents the two particle irreducible quark-antiquark potential, and the round 
blob is the full T-matrix. The lines correspond to constituent quarks ( thick lines 
in Fig. (0)). 



3.5. The pion as quark-antiquark bound state 

In the spontaneously broken phase, the constituent quark mass, M 7^ 0. 
As a consequence of Goldstone's theorem there must exist massless pseu- 
doscalar particles, which one identifies with the pions. This can actually 
be obtained by solving the Bethe-Salpeter equation in the pseudoscalar- 
isovector channel for quark-antiquark scattering (see Fig. ^) and checking 
that there is a bound state pole which becomes massless as the current 
quark mass goes to zero, m ^ 0. In the lowest order approximation, where 
the kernel is approximated by a constant, the Bethe-Salpeter amplitude for 
vr"*" becomes 

Tap;^s = (75r+),^ (75r+)^5 t{P) (39) 
where t(P) is a number depending only on the total momentum, given by 

1 1 



75- 



M 



t{P) (40) 



In writing the polarization operator, a particular choice of momentum rout- 
ing has been made, but in a gauge invariant regularization such as PV, one 
may safely shift the integration variable by any amount of the external mo- 
mentum P. Euclidean 0(4) cut-offs do not preserve this property. Using 
the PV regulators one gets 

t(p)-i = - 8NJ + 4AfcP^F(/) (41) 

where / is defined by Eq. (|29| ) and the one- loop function, F{p'^), is given by 

d'^k 1 1 



F{p') = (-i) 



(27r)4 A;2 - M"^ {k - pY - 



(42) 
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Fig. 3. Pion pole contribution to the Bethe-Salpeter quark-antiquark scattering am- 
plitude. The small solid blob represents the Bethe-Salpeter pion quark-antiquark 
amplitude. Solid lines are constituent quarks. . 



Using the Feynman trick for the two propagators, this function can be 
rewritten as 



dxF{p'^, x) 



(43) 



where 



F{p^,x) = (-i) 



1 



(47r)2 



(27r)4 [k2 - M2 + x(l - a;)p2 + ie]' 
- Ci log + Af - x{l - x)p" 



(44) 



In the chiral limit, m = 0, the gap equation ( |35| ) can be used, and the first 
terms in the r.h.s. cancel, producing a massless quark-antiquark bound 
state. However, for p"^ > AM"^, the function F{p'^), and hence the t-matrix, 
develops an imaginary part indicating a lack of confinement. For <C 4M2 
we get a deeply bound state, and thus hope confinement not to be essential. 
For M = 300MeV one has ml/{4M^) ~ 0.06. 



For a finite quark mass we have a pole at p 



ml 



m 
GM 



4NrmlF(ml: 



(45) 



The pion coupling of a composite pion state to a quark-antiquark pair is 
given by the residue of the t-matrix at the pion pole (see Fig. |^ ) , 



9: 



dt-^{p^ 



dp^ 



4N, 



d[p^F{p'^)] 



dp^ 



So that the composite canonically quantized pion propagator is 

Kp") 1 



Dip' 



lirqq 



P^ 



mt 



{p^ 



(46) 



(47) 
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Fig. 4. Pion (dashed line) weak decay into an axial current (wavy line). The solid 
blob represents the Bethe-Salpeter pion quark- antiquark amplitude. Solid lines are 
constituent quarks. . 



3.6. Pion weak decay 

Using the bound state solution of the Bethe-Salpeter equation for the 
pion, one can compute the pion weak decay defined as 

{0\jfm7r''ip))=iUp'^6''' (48) 

yielding (see Fig. |), 

U = 4N,Mg^ggF{ml) (49) 

and thus 

which, in the chiral hmit yields the Goldberger-Treimans relation at the 



quark level, g-Kqqfn = M. Combining Eqs. (pq), (45) and ( p9[ ) one gets 



-m{gg)=mlfJ^L^] (51) 

which at lowest order in the chiral expansion becomes the well known Gell- 
Mann-Oakes-Renner relation. 

Eqs. (|49|) and (|^) are usually employed to fix the parameters in the 
broken phase, taking for /^r = 92MeV and m-,^ = 139.6MeV. For a given 
constituent quark mass, M, one can compute the cut-off A. The gap equa- 
tion allows to determine the coupling constant G. For instance, using the PV 



regularization given by Eq. (|3|), for M = 280MeV one obtains A = 871MeV 
and {qq) = — (290MeV)^. From this number one may deduce that, since it 
resembles the numerical value of the condensate at the scale /i = IGeV, 
that the model scale is around IGeV. Of course, the numerical values may 
change in a different version of the PV scheme (like e.g. with Ai 7^ A2, or 
adding more terms in the sum) . Fortunately, as we will see below, the results 
for some high energy properties are insensitive to the choice of parameters, 
up to unimportant chiral corrections. 
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Fig. 5. Electromagnetic pion form factor in the NJL model. The closed loops are 
Pauli-Villars regularized. Wavy lines arc photons, solid lines are constituent quarks 
and dashed lines composite pions. 

3. 7. Pion wave function 

In principle, the Bethe-Salpeter pion-quark-antiquark wave function is 
given by (we take vr"*") 

^^^^^ = f/2 +\ - M (-^-'^'^^^^^) r/2-\-M ^'^^ 
but PV regularization has to be understood because, at least formally, 

UP, = (-i) J ^Tr (^XP(A:)^7^75) (53) 

and the integral is logarithmically divergent if Eq. (^) is taken literally. 
This problem is typical in chiral quark models and has been the cause of 
much confusion. Regularization is only easily applied at the level of closed 
quark lines 0. We will see more of this later, but it has always been a weak 
point in quark- loop calculations. One obvious choice to regularize open lines 
is to enforce consistency with the closed lines, although at first sight this 
procedure looks a bit arbitrary. Another possibility to open a regularized 
quark line, but starting from a closed quark line, is to consider a physical 
process involving photons and pions in the high energy limit and extract 
the leading power behavior. Prominent examples which will be analyzed 
below are deep inelastic scattering from which the pion structure functions 
can be deduced in the Bjorken limit or the pion transition 7* 7r''7-form 
factor in the limit of high photon virtualities and fixed photon assymetry 
from which the pion distribution amplitude may be derived. 

3.8. Pion Electromagnetic Form Factor 

The pion electromagnetic is probably the simplest case where one can 
illustrate some of the points we want to make. For a charged pion vr"*" = ud 

^ This is actually a good reason to use bosonization schemes; there the concept of an 
open quark line does not appear . 
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the electromagnetic form factor is defined as 

{7r+{p')\j;^m^^ip)) = err(p',P) = ip^+p'neFriq') (54) 

and can be computed using Fig. ^. For on-shell pions and in the chiral hmit 
the result is rather simple 

Jit 

where F{p'^) is given by Eq. (p2|). The pion form factor is obviously properly 
normalized, i^^™(0) = 1, as it should be due to the gauge invariance of the 
Pauli-Villars regularization. The mean-square pion radius reads 

a result explicitly depending on the regularization. Numerically one has a 
reduction of about 25% for M = SOOMeV due to the finite cut-off effect, 
((r2)|™)V2 = (o.58fm) x 0.75 = 0.50fm, to be compared with the experi- 
mental number, 0.66fm. Two points are worth stressing here. Firstly, it is 
clear that the finite cut-off corrections go in the wrong direction. On the 
other hand, one should not expect a perfect agreement with the experimen- 
tal number, since one expects a sizeable contribution from the pion loops 
{l/Nc corrections in the model), as suggested by ChPT There is a way 
one can get rid of the cut-off corrections; one can simply split the charge 
contribution at q'^ = 0, which should be one, and compute the difference 
using a unregularized quark loop. This is equivalent to write down a once 
subtracted dispersion relation for F^^{q'^). For the unregularized quark 
loop only one has the asymptotic behavior for — oo, 

Friq')^\ogi-qyM^), (57) 

without power corrections. For the PV regularized NJL model, at high 
Euclidean momentum one has instead a cut-off independent relation, which 
is a pure power correction ^, 

q'Friq') - ~ 0.34GeV2 (58) 



for M = SOOMeV and the PV method, Eq. (|3|). Up to the measured 
momentum transfers, — 2GeV^ < q^ < — 6GeV^ this value is very similar to 

^ For a Euclidean 0(4) cut-off one gets a logarithmic behavior log((5^/A^)/Q^ |H. 
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the experimental average number 0.38 it 0.04 |^^. However, it is in principle 
not clear if this number should directly be compared to experiment, without 
taking into account QCD radiative corrections. Another point of view is to 
reject the validity of the model to such high energies, matching instead to 
the known QCD result at LO g (see also Ref. ||], 

q^Fr{q^) - lQT^a{-q^)fl (59) 

which for cp' = — 4GeV^ yields, 0.13, a too small value as compared to the 
experimental one 0. We may identify both coefficients at the model scale, 
//o, yielding the result 

In this relation we have stressed the fact that also the constituent quark 
mass is fixed at the model scale, //q- Actually, we may use relations ( |6H ) 
and ([l5|) to get, 

M^o = 300 ± 80MeV (62) 

Although this value should be considered a crude estimate it provides a rea- 
sonable value for the constituent quark mass at the model scale. Conversely, 
if we assume M = SOOMeV then a(//o)/(2vr) ~ 0.5 it 0.2, a compatible value 
with the momentum fraction estimate. 



4. Pion Distribution Functions 

Deep Inelastic Scattering 

Inclusive lepton-hadron scattering is described in terms of the hadronic 
tensor, W^u{p,q), and can be obtained from the imaginary part of the for- 
ward Compton amplitude for virtual photons as follows [11 1 



zvr 



Wi{q\p-q) [-g^,u + 



+ 



W2{q^,p 



mt, 



Actually, this expression can be obtained from 



p-q 



^ ipT,{x, Q) 
6a;(l — x) 



(63) 



Pu 



P-. 



(60) 



with ipTr{x,^) the pion distribution amplitude at the scale ^, when the asymptotic 
wave function, 99^(2;, cx)) — &x{l — x) is substituted H (see also below). 
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Fig. 6. The Compton amplitude in the Bjorken hmit and its relation with the quark- 
target scattering in the parton model. The high and low energy contributions are 
also explicitly displayed. 



where 



(64) 



Here, rrip = is the mass of the pseudoscalar meson and q is the mo- 
mentum of the virtual photon. In the Bjorken limit, one obtains [11| 
ix = -qy{2p-q)), 



Wi{x,Q^) ^ Fi{x) 
W2ix,Q^)^F2{x) 



(65) 
(66) 



where the Callan-Gross relation F2{x) = 2xFi{x) is fulfilled, as a conse- 
quence of the spin 1/2 nature of quarks. The structure function is given 

by 



(67) 



i=uA 



where = 2/3 and = —1/3 are the quark charges and qi[x) and qi{x) 
are the momentum fraction distribution of the different quark species. 

Introducing light-cone (LC) coordinates x = (x+,x~,x^) with x^ = 
x^ lb the analysis of Ref. [|ll[ yields 



Fi{x) 



4tt 



dk <fk± 
(27r)3 



tr 



(68) 



where Q = diag(e„, e^) is the charge operator and the forward quark-target 
scattering amplitude is defined 



X{P, k) 



dHe'^Hp\T{qm{mP) 



(69) 
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xip^k) corresponds to the unamputated vertex. Eg. ([6^) holds under the 
assumption of scaling and finiteness in the Bjorken limit. 

Logarithmic scaling violations due to perturbative QCD radiative correc- 
tions , relate in a linear fashion the leading twist contribution to structure 
functions at a given reference scale, /io, to the scale of interest, /i, (in DIS 
one takes Q'^ = fj?) 

Fl{x,^x) = u{^l,^lo)Fl{x,^lo). (70) 

where Eq. (^) is satisfied for the linear matrix integral operator U. Asymp- 
totically, QCD predicts Q 

F2(x,Q2) = 2xFi(x,Q2) ^__^^d__s(^x) (71) 

Actually, the pion structure functions cannot be measured directly since 
there are no pion targets. Instead, the parton structure of the pion has 
been analyzed [l^ and a simple parameterization at = 4GeV^ has been 
given. The valence quark distributions extracted in Ref. |jl^ from Drell-Yan 
experiments [^] seem well determined, whereas the gluon distribution as 
obtained from the analysis of prompt-photon emission data [37| is less well 
determined. A recent analysis p^] of the ZEUS di-jet data seems to favor 



the gluon distributions of Ref. [IS]. This is why, in determining the low- 
energy scale of our model, we use the valence momentum fraction of 47% 
found in Ref. p| at AGeV"^ instead of 40% found Ref. pi ^] using some 



additional assumptions |56]. 



A result found long ago and usually adopted in DIS calculations are 
the so-called counting rules. They establish that there is a relation between 
the asymptotic behavior of the electromagnetic pion form factor for large 
Euclidean momenta and the structure function behavior as x ^ 1, namely 



Friq') ~ l/{q'r Flix) ~ (1 " 2^)'"-' (72) 

In QCD F^"^{q'^) has a logarithmic radiative correction, Eq. ( |59D so we may 
take n ~ 1, which implies F{x) ~ 1 — x. Although it is not clear at what 
scale the counting rules are valid, recent investigations based on quark- 
hadron duality confirm the exponent n ~ 1 for F{x) at /U ~ 2GeV 



4-2. Spectator model and quark loop calculations 

The spectator model is a very simple model where one can illustrate 
the source of the problems in this kind of calculations and the relation of 
gauge invariance and normalization. Actually, in the case of the pion such 
a model corresponds to a unregularized NJL calculation. The spectator 
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model consists of replacing the full sum of intermediate states in the soft 
piece of the Compton amplitude by a spectator particle, which for the pion 
has quark quantum numbers All is required to do the calculation is to 
know the vrgg-coupling. We take a pseudoscalar coupling as follows 

Kqq = -l^r'^g^m (73) 

At the level of the Compton scattering amplitude it corresponds to use the 
well-known hand-bag diagram. Straightforward calculation of the unregu- 
larized structure function yields in the Bjorken limit 



- xil - x)ml 



log 



M2-x(l-x)m2 J ^ ' 



which is a scaling violation quite different from those found in QCD. Actu- 
ally, we get F2(x,oo) = 00 instead of Eq. ([Til) . Scaling can be restored by 
attaching a form factor to the nqq vertex or putting a 0(4) cut-off in the 
loop. No solution is really free of problems, since either gauge invariance is 
violated or extra singularities are introduced. 

At the level of the quark-target scattering amplitude the spectator pic- 
ture corresponds to consider an intermediate constituent quark state. The 
result for a pseudoscalar coupling is 

Fiix) = 7T(e„ + ea) / ,o^^2 r,.2 , .^2 33212 (^5) 



2' " 47r 7 (27r)2 [A:2 + M2 - x(l - x)m2]2 

However, we see that the integral in k± is actually logarithmically diver- 
gent. This only reflects the fact that for the unregularized quark loop the 
separation between soft and hard processes involved in the Bjorken limit 
takes place at /c^ ~ Q2_ Obviously a form factor for the vertex or a trans- 
verse cut-off, \k±\ < A (both independent), may be introduced yielding 
a finite result. 

Thus, we learn in this particular model that the connection between 
the Compton amplitude and the quark-target scattering formulas is only 
justified provided everything is finite. So, some high energy suppresion 
should be applied, with the help of a regularization. The problem is that 
we expect to do it in a gauge invariant manner, since after all we are dealing 
with an electromagnetic process. 

4-3. PDF Calculation from the Compton amplitude 

In the PV regularized NJL model, the gauge invariant forward virtual 
Compton scattering amplitude is given the sum of the hand-bag diagram. 
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Fig. 7. Spectator model. The sum over intermediate states in the soft hadronic 
contribution is replaced by a single quark state. 



the crossed contribution and the 77r — > vr* ^ tt^ involving the ofF-shell 
electromagnetic pion form factor (see Fig. (|8|)). Actually, it can be shown 
that in the Bjorken limit only the hand-bag contribution survives, since the 
crossed and pion off-shell contributions are higher twist. The result has 
been obtained in Refs. |p5| , p^ j and turns out to coincide with applying PV 
regularization, Eq.([3l]), to Eq. ( |7^). The result is 



1 



d 



where the function F(x,p^) is given by Eq. 
one gets in the chiral limit 



(76) 



. Taking 7r+ for definiteness. 



d^+{x) = 1 , < X < 1 . 



(77) 



This means the following result for the valence-quark distribution function 



V{x) 



< X < 1 



(78) 



As we see the result is independent of the PV regulators and has proper 
support, normalization and fulfills the momentum sum rule. Moreover, 
Eq.(|78|) is consistent with chiral symmetry, in the sense that it has been 
obtained taking explicitly into account the boson Goldstone nature of the 
pion. The result has also been obtained within the NJL model imposing a 



transverse cut-off [42| upon the quark-target amplitude ([75|). 

Eq.(|7D disa grees with other quark loop calculations. If Eq.(|64[) is used 
with a four-dimensional cut-off |^3[ or Eq. ( 65 ) is used with Lepage-Brodsky 
regularization |42] in the NJL model, different shapes for the quark dis- 
tributions are obtained. The null-plane [44|, NJL model ||4^, spectator 
model [41 1 and the recent quark loop calculations also produce different 
results. In all cases, the use of momentum dependent form factors or non- 
gauge invariant regularization make the connection between Eq.(|6^) and 
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2^ 




Fig. 8. Gauge invariant Compton amplitude in the NJL model. The closed loops are 
Pauli-Villars regularized. Wavy lines are photons, solid lines are constituent quarks and 
dashed lines composite pions. 



Eg .(1681) doubtful and, furthermore, spoil normalization. The results based 
on a quark loop with momentum dependent quark masses |51, ^] seem to 
produce a non-constant distribution. 



4-4- PDF Calculation in Light-Cone Coordinates 

Perhaps the simplest way to obtain PDF for the pion in the NJL model is 
by using the pion Bethe-Salpeter propagator. The probability distribution 
is defined through the pion expectation value of the of number operator 
corresponding to a quark of a given species. Taking into account the proper 
normalization factors one gets at the pion pole the following expression 
for the probability distribution of finding a valence quark with momentum 
fraction x and transverse momentum k± 
1, 



-V^{x,kjL) = / dk-D{p,k) 



fc+=m,ra;,p^=m,r 



(79) 



where the Bethe-Salpeter normalization in the Pauli-Villars regularization 
scheme is 

d r 9^ 1 1 I 

(80) 



D{p,k)=m,igl 



Ci- 



k^-M^-Af + ie {k - pf - M2 



Introducing light-cone (LC) variables 

— — — 

Employing the definition ( |79|) one finds 



d^k 



1 



-dk+dk-d^k± 



(81) 



^V^{x,k±) 



2iN,gl 



qq 



dk~^dk 



(27r)4 x{l-x) 
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5 {k~^ — xp'^ 



1 



k- 



fc2+M2+A2+iO+ 



fc2 +M2+A2+iO+ 



where the location of poles in the k" variable has been explicitly displayed. 
For X > 1 or X < both poles are above and below the real axis respectively, 
and hence the integral vanishes in either case. For < x < 1 the integral 
yields 



^V^{x,k_i 



nqq 



d 



16tt^ dml 



' k'j_ + m? + A? — m^x{\ — x) 



.(82) 



Notice that, due to the Pauli-Villars subtractions, we have the asymptotic 
behavior in the transverse momentum k \ , 



-V;(x, A;_l) 



167r3 k^, 



(83) 



This guarantees the convergence of the k± integral without introducing a 
transverse cut-off. Thus, the Pauli-Villars regulators automatically provide 
a (gauge- invariant) form of a transverse cut-off. Integrating the transverse 
momentum we get the PDF 

1 



:K(x) 



1 



d k±-VT,{x,k±) 



d 



mlF{ml,x) 



(84) 
(85) 



which is the result found in Ref. |^ (see also Eq. (|7|) above). In the 
chiral limit, m-j^ = 0, one can use the Goldberger-Treiman relation for the 
constituent quarks, gnggfn = Then = 4A''cM^F(0), which gives the 
very simple formulas 



^V-^ix,kj_) 



2^-(^) 



1. 



1 



(86) 
(87) 



In the chiral limit V^(x,A;j_) becomes trivially factorizable, since it is inde- 
pendent of x. A remarkable feature is that the last relation, Eq. (^), is 
independent of the PV regulators. 

The light-cone interpretation has been pursued (see e.g. [44| and refer- 
ences therein) and more recently |69, 42] within a LC quantization. In these 
cases transverse cut-off's were introduced, a posteriori. As we have shown 
above this is not necessary in the PV regularization. 
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4-5. QCD evolution 

Due to radiative corrections, parton distribution functions evolve loga- 
rithmically with scale through the DGLAP equations |^. Non-singlet (or 
valence) distribution functions are easily evolved in terms of their Mellin 
moments. For vr"'" = ud one defines 

V{x) = u{x) - u{x) + d{x) - d{x) (88) 

and takes (we assume chiral limit for simplicity) ^ 

^K(a;,/xo) = l (89) 
At leading order, evolution reads 

Vnili)^ dxx^V{x,,i)=(-p-] / dxx"y(x,^o) (90) 

Jo \a[iJ,o)J Jo 



where the anomalous dimension is 

7n = -^Cf 



2 1 

n + 1 n + 2 41^; 



(91) 



with Cp = 4/3. Taking n as a complex number, which also requires an 
analytical continuation of both V^(/Uo) and 7^^, Eq. (|9^) can be inverted 
using 

pc+iao Aj^ 

y{^,i^)= ;r^^-"K(/^) (92) 



27ri 

where c has to be chosen as to leave all the singularities on the left hand side 



of the contour. The result for the pion valence PDF found in Ref. |25| at LO 
and |2|1 at NLO is displayed in Fig. 1^ ai ^'^ = = 4GeV2 compared with 
phenomenological analysis for the pion As one can see, the agreement 
is quite impressive. Despite the fact that a(^o)/(2vr) ~ 0.3 the differences 
between LO and NLO turn out to be small. 

It is interesting to analyze these results on the light of the counting rules, 
Eq. (^), of Ref. [39|, as was already done in Ref. p^. Using the inverse 



For non experts it may sound unnatural to take the high energy, 00, hmit in a 

model and use it as initial condition for QCD evolution at a low scale /io. The point is 
that if the quark model scales, the asymptotic behavior can be separated according to 
increasing power corrections in 1/Q^ in a twist expansion. The anomalous dimensions 
relevant for QCD evolution of structure functions are only known for the lowest orders 
of such an expansion, but not for the full structure function. Using Eg. (^) as initial 
condition automatically complies with the asymptotic result of Eq. (\ 
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Fig. 9. Valence distributions in the pion at Q4 = 4GeV^ at LO p| and NLO M 
compared with phenomenological analysis for the pion SMRS92 |13[, GRV99 pSj 
We take (xV)^ = 0.47 at = AGeY^ M. 



Mellin formula, Eq. (92), it can be shown that, if V{x,iJ,o) 
then 

y(x,/x) ^ C(l -X) —'"TO 5 



C{1 



\N 



X 1 (93) 
The value of the additional contribution to the exponent is a weakly de- 



pending function for > IGeV. For the pion we find, using Eq. (15) 



y(a;,2GeV)^2(l 



il.liO.l 



1 



(94) 



which is consistent with the counting rules and the phenomenological anal- 
ysis 1 13 1, after evolution. A more detailed account of the sea and gluon 
distribution functions at LO and NLO and also for K and r] mesons can be 
found in Ref. 



5. Pion distribution amplitude 

5.1. Pion Transition Form Factor 
The matrix element for the 7* 7r*^7 transition form factor is defined 

r(^o_^^^(p,9i,g2) = W/39l92-^^77(^''9l'92) (95) 
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where -^7777(^,^1,^2) is the vr'^ electromagnetic transition form factor Q. 
Two QCD results apply for this transition form factor. In the chiral limit 
and for on-shell pions, = 0, and photons, Qi = q2 = 0, corresponding to 
the t:^ — > 27 decay one gets the normalization condition 



1 



F,^^(0,0,0) = (96) 



which is the standard result expected from the chiral anomaly. The defi- 
nition of the leading twist contribution to the pion distribution amplitude 
(PDA), ip.^{x), from the transition form factor is 



lim Q2 F„yy{p,qi,q2] 



,2 = _(i+^Q2^,2^_(i_^Q2 



= r dx ^-(^^ (97) 

Nc Jo 1 - w(2x - 1) ^ ' 

where to = {ql — ^2)/(9i + ^2) photon asymmetry, — 1 < lj < +1, 

which is kept finite when taking the limit ^ 00. It can be shown that 
in the parton model the PDA can be computed by 

iPn{x)= J (fk^^^{x,k±), (98) 

where the pion light-cone wave function (pseudoscalar component) is defined 
as the low-energy matrix element ||4^, ^] 



{n+{p)\u{r,Uh^l5dmO). (99) 
where = m-j^ and p± = 0± . An important relation found in Ref. Q reads 

U C dx'^^{x, 0^) = F^77(0, 0, 0) (100) 
Jo 

Radiative logarithmic corrections to the pion distribution amplitude (PDA) 
can be easily implemented through the QCD evolution equations Q, which 
yield for Q"^ ^ co the asymptotic wave function of the form 

ip-n{x^(X)) = Qx{l — x). (101) 

Moreover, the pion transition form factor has been recently measured by 
the CLEO collaboration [^] and a theoretical analysis of PDA based on 
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these data and light-cone sum rules has been undertaken , showing that 
at Q = 2.4 GeV PDA is neither asymptotic, nor possesses the structure 
proposed in early works |]60|| . 

The pion distribution amplitude has been evaluated with QCD sum 
rules |^l|], in standard [ p^] (only the second ^-moment) and transverse lat- 
tice approaches [63, 54, |65f|, and in chiral quark models Il66|, |67l 
71 



70 



2|, 74|. In chiral quark models the results are not always compatible 
to each other, and even their interpretation has not always been the same. 
While in same cases there are problems with chiral symmetry and proper 
normalization [36, 67, in other cases [68, [HT 



70, 71 



, |7^, |7J[ it is not 

clear how to associate the scale at which the model is defined, necessary to 
define the starting point for the QCD evolution. The fact that several cal- 
culations [^, ^ |7^, 0, |73| produce a PDA strongly resembling the 
asymptotic form suggests that their working scale is already large, and the 
subsequent QCD evolution becomes unnecessary, or numerically insignifi- 
cant. This also tacitly assumes that these models already incorporate the 
QCD radiative corrections. 



5.2. PDA calculation for the unregularized quark loop 

Since 7* — > -k^^ is an abnormal parity process, the standard procedure 
in the NJL model is not to regularize it because this is the only way to 
preserve the anomaly (See also Ref. p5| , [7^])). Straightforward calculation 
yields the result 



1 



1 



1 



i J (2^)4 A;2 - M2 {k 



Qi'j 



M2 (k 



(12) 



^1 



^2) 



In the chiral limit g-wqqf-w 

q2 = 0, corresponding to the vr'^ 



1! 



M and for on-shell pions, = 0, and photons, 
27 decay one gets 



7777 



(0,0,0) 



8M2 
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d^k 



1 



1 



(27r)4 (A;2 - NPf Att'^U 



(103) 



which agrees with the chiral anomaly expectations, Eq. (|96|). This result is 
the main motivation for not introducing an explicit cut-off^ in the abnormal 
parity processes. In order to compute the the high momentum behavior, 
we use the Feynman trick in the two propagators containing qi and ^2, and 
shift the integrating variable. For on-shell massless ir^, p'^ = 0, we get 



F * . 



8M2 

1 



d'^k 



1 



(27r)4 (fc2 



M2)2 

1 



dx 

{k + xqi + (1 



x)q2)^ - M2 



(104) 
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Fig. 10. Transition form factor in tiie NJL model. The closed loop is Pauli-Villars 
regularized. Wavy lines are photons, solid lines are constituent quarks and dashed lines 
composite pions. 



The — > CO limit can be undertaken yielding (see Eq. (| 



(47r)2 



log 



- x{l - x)ml 



(105) 



where a finite pion mass has been reinserted. Thus, the prescription not to 
regularize the abnormal parity vertex does not agree with the factorization 
result found in a parton model approach. Instead, there appear scaling vi- 
olations, which do not correspond to those expected from QCD evolution, 



contradicting the asymptotic result, Eq. (101). This problem is of a simi- 
lar nature as the one found in the discussion of the spectator model after 



Eq. (74), and similar remarks concerning gauge invariance apply here. 



5. 3. PDA calculation from the transition form factor 

To achieve factorization, it is necessary to effectively cut-off the trans- 
verse momentum, i.e., we have to be able to make larger than any scale 
in the loop. This is not consistent with integrating in k up to infinity^. If 
we introduce PV regulators [|, we get 



16M2 ^ 1 r d^k 1 

'J^^'^'lJ (2^)4 (fc2-M2-A2)2 



1 



X / dx -, r 

/o 1 - uj{2x - 1) 



(106) 



^ If we formally take the limit we get an expression looking like Eq. ( p7| ) with, = 1 

but with the unregularized form of ( see Eq. (^). To achieve Eq.(^) with the 
regularized definition of some regularization must be introduced. 

^ This regularization guarantees that the momentum routing is irrelevant. The Eu- 
clidean 0(4) cut-off of Ref. (68[ Q requires a very special momentum routing after 
regularization; a different momentum routing with the same cut-off would produce 
finite cut-off contributions to their results. 
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which, using Eq. (49) and Eq. ( |97| ) corresponds to take (after inclusion of a 
finite m.^ ) 

if^ix) = 4N,M^F{mlx), (107) 

Jtt 



where F{p^,x) is given in Eq. (^^. ipT^{x) depends only on x and is properly 
normalized (see Eq. (^)). In the chiral limit one gets 

if^{x) = l 0<x<l. (108) 

independently on the regulators. However, if we PV regularize the vertex 
there is a violation of the anomaly, since 

1 

F.,,(0,0,0)|,,g = ^^c.^^^^ (109) 

For typical values of the parameters one finds a 25% reduction in the am- 
plitude which means a 40% reduction in the width for vr'^ — > 77. Such a 



result was also found in Ref. |35| using a Euclidean 0(4) cut-off. Thus, in 
the present framework we have to choose between preserving the anomaly 
and not obtaining factorization or vice-versa, i.e. violating the anomaly and 
reproducing factorization. 

There is a subtlety in the previous reasoning, because the prescription of 
not regularizing the abnormal parity contribution to the action really means 
that the result is conditionally convergent and can be given unambiguously 
if one insists on maintaining vector gauge invariance [^] . A practical way to 
implement this fact is to introduce a gauge invariant regulator and remove 
it at the end of the calculation. The triangle graph is linearly divergent, and 
thus a regularization must be introduced. If one insists on preserving vector 
gauge invariance the regulator must preserve that symmetry, but then the 
axial current is not conserved generating the standard chiral anomaly. The 
obvious question arises whether the limit — > cxd must be taken before 
or after removing the cut-off by sending it to infinity. If one takes the 
sequence > — > 00, the results of Ref. produce a constant PDA, 
in agreement with our low energy calculation. 

5.4- PDA calculation in light- cone coordinates 
The pion distribution amplitude is defined through 

r rlh+ dh^ d'^h 

Up^^^{x) = J ^(^^)4 ^ (^^ - ^P'') [XP(A:)7^75] (110) 
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where the PV regularized form the pion BS wave function, Eq. (|52[), is 
understood. Formally, in momentum space, Eq. (|9^ ) corresponds to inte- 
gration over the quark momenta in the loop integral used in the evaluation 
of /tt Eq. (^9|), but with = p^x = itIt^x and k± fixed. Thus, with the PV 
method and after working out the Dirac traces, we have to compute 

1 1 



A;2 +M2+A2+iO+ , A:2 j^M'2+k'^.+iO+ 

ryy-t 2. 3 

A; / / b'jr /-I \ Irj 

" mTr(l—x) m-^x 

where, again, the location of poles in the k~ variable has been explicitly 
displayed. This integral coincides with that found when computing the 
PDF (see Eq. (p2[)). Evaluating the k~ integral gives the pion LC wave 
function in the NJL model with the PV regularization: 

ANcMg^rgq ^ 1 

"^"^"^'^^ = 167r3/, ^''^■/ti+A2 + M2-x(l-x)m2- ^^^^^ 

For m,r 7^ it is non-factorizable in the k± and x variables. As a consequence 
of the PV condition with two subtractions one has, for large k±, 

^.{xM) - ^^^^37^ (113) 

which gives a finite normalization and a finite second transverse moment, 

{ki) = J (fk± J\x^^{x,k^)kl. (114) 

Integrating with respect to k± reproduces the pion distribution amplitude of 
Eq. ( [1071 ). In the chiral limit, = 0, one can use the Goldberger-Treiman 
relation for the constituent quarks, g-nqqf-K = yielding in addition to 
Eq. ( |108| ) the very simple formulas 

^7r(a;,fc±) = ]^V^{x,ki_), (115) 



.2 , _ M{qq) 



(116) 



where ^-(a^, k±) is given by Eq. ([Sq). Note that these relations are indepen- 
dent of the PV regulators. 
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One aspect of regularization should be analyzed here. According to 
previous studies @], the value of the LC wave function in the chiral limit is 
fixed at k± = 0_l by the chiral anomaly, Eq. ( |100| ). In the PV regularized 
NJL model this is not the case, since 



/■I AM m"^ 1 

Ul <ix>l.,(.,OJ = F(0,0,0) = -^^c.-^^ (117) 



The first two terms in this equation indicate the consistency of our calcula- 
tions between the transition form factor and the low energy matrix element, 
but confirms the anomaly violation we have referred to above. 

This is a clear deficiency of the NJL model and its regularization pro- 
cedure, and we do not know of any convincing way of avoiding this prob- 
lem in this model ^ Nevertheless, it can be shown that it is possible to 
write down a chiral quark model based on the concept of spectral reg- 
ularization where the proper anomaly is reproduced and the results 
fn{x,lJ-o) = VTr{x, fj,o)/2 = 1 still hold 



5.5. QCD evolution 

The leading-twist PDA requires, the inclusion of radiative logarithmic 
corrections through the QCD evolution For the pion distribution ampli- 
tude this is done in terms of the Gegenbauer polynomials, by interpreting 
our low-energy model result as the initial condition. In the chiral limit, 

l.p■n{x,^lo) = 1. (118) 

Then, the LO-evolved distribution amplitude reads Q 

oo ' 

(^,(x,/x) = 6x(l-x)^ C73/2(2x-l)a„(^), (119) 

ra=0 

where the prime indicates summation over even values of n only. The matrix 
elements, ttnifJ-), are the Gegenbauer moments given by 

2n + 3 / a(^y"'/(2/5°) 



""^^^ - 3(n + l)(n + 2) U(/io)j 

dxC^/^i2x-l)ip^ix,fio), (120) 







In the Euclidean version of the model the accepted regularization prescription is to 
regularize the real part of the action (normal parity processes') and not to regularize 
the imaginary part of the action (abnormal parity processes) ItsI . Such a prescription 
agre es wi th the anomaly but does not agree with factorization and produces instead 



Eq. (105) 
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X 

Fig. 11. The pion distribution amplitude in the chiral hmit evolved to the scale 
Q"^ = (2.4GeV)^. The two values for the evolution ratio r ~ oi{Q)/a{Qo) reflect 
the uncertainties in the values of Ref. based on an analysis of the CLEO 
data. We also show the unevolved PDA, (Ptt{x, Qq) = 1, and the asymptotic PDA, 
(fTr{x,oo) — 6x{l — x). We use, as suggested by the analysis of the DIS, a{Qo) — 
2.14 of Ref. H, which corresponds to Qq = 313 MeV. 



with Cn^'^ denoting the Gegenbauer polynomials, and > is given by 
Eq. (|l]). From Eq. ( pSj) one gets immediately 

[\xC^/\2x-l)^^ix,Qo) = l. (121) 
Thus, for a given value of /U we may predict PDA. To determine the initial 



scale fiQ, or, equivalently, the evolution ratio r = a{Q)/a{Qo) Ref. ||2^ uses 
the result of Ref. ||, where it is found a2(2.4GeV) = 0.12 ± 0.03. Using 
this input, one gets 

fio = SlStfo^eY (122) 

Within uncertainties, this result is compatible with the values obtained from 
the momentum fraction analysis, Eq. ( p!5| ) and the pion electromagnetic form 
factor, Eq. (pTj). The result for the pion distribution amplitude obtained in 



Ref. |29| is shown in Fig. 11 for Q = 2.4GeV and reflecting the uncertainties 
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from Ref. |5£]. After evolution the results closely resemble those found in 
the transverse lattice 34, 35 1. The analysis of the end point behavior 
yields, after evolution, the estimate 



ACf , a(/i) 



In 



/5o a(/^o) 



+ 1 



12.5x. 



(123) 



for X — > 0. Here C,{s) = is the Riemann zeta function and the 

numerical value corresponds to take ^ = 2.4GeV. The value obtained in 
Ref. [29 1 for the second ^-moment (^ = 2a; — 1) is 



(<f2) = / dx(p^(x,g = 2.4GeV)(2x-l)2 
Jo 

= 0.040 ±0.005, (124) 

to be compared with (^^) = 0.06 it 0.02 obtained in the standard lattice 
QCD for Q = l/a = 2.6 ± O.lGeV ||. 

Based on the identities v9jr(x,//o) = K-(2;,Aio)/2 = 1 between the PDF 
and the PDA in the chiral limit at the quark model scale a remarkable 
integral relation between (p-,^{x,fi) and V-n{x^^)/2 has recently been derived 
using LO evolution equations l2r 

9^77(3;, /i) 
6a;(l — x) 

where the explicit expression for the scale independent kernel, K(x,y) is 
given in Ref. [29|. Using this equation one can regard Fig. 11 as a predic- 
tion of the PDA in terms of the PDF parameterizations of Ref. ||l^ (see 
Fig.D. Using Eq. ( pSj) one gets the following estimate for the leading twist 
contribution to the pion form factor at LO 

dx^4^ = 1.25 ± 0.10 (126) 



1 



(125) 



2A 



Twist- 2 



6x(l — x) 



The experimental value obtained in CLEO ||5^ for the full form factor is 
Q'^F^*^^^{Q)/2f^ = 0.83 ±0.12 at = 5.8GeV^ Taking into account that 
we have not included neither NLO effects nor an estimate of higher twist 
contributions, the result is rather encouraging. Finally, using this estimate 
for ipTj{x,Q) in the electromagnetic form factor, Eq. (^), one deduces an 
enhacement 



Q^Fn-Q') 



Twist -2 



0.20 ± 0.03 



Q2 = 5.8GeV2 (127) 



which also brings the number closer to experimental number [^], although 
it only accounts for half of its value ~ 0.38 ± 0.04GeV^, perhaps due to 
the lack of NLO or higher twists. At present this point seems not to be 
understood and deserves further investigation. 
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6. Conclusions 

Chiral quark models incorporate two essential features of QCD at low 
energies: spontaneous chiral symmetry breaking and quark degrees of free- 
dom. In addition, it can be made compatible with large Nc counting rules, 
and calculations in the leading order approximation have been undertaken 
in the past. In these lectures we have focused and reviewed particular appli- 
cations in a prototype chiral quark model, the Nambu-Jona-Lasinio model, 
but many results extend trivially to other models. We have also restricted 
to the pion because we do not expect a better theoretical understanding 
of a hadron, including the fact that confinement is hoped not to play an 
essential role. 

A very important issue concerning the treatment of low energy models 
is the existence of a high momentum suppression in the interaction. In the 
NJL model this is done via a regularization method, which has to comply 
with several properties, like gauge invariance and scaling in the high energy 
limit, i.e. the absence of spurious logarithmic corrections. This allows to 
clearly identify power corrections of pion observables. 

These models make sense at a given low renormalization scale /iq. Any 
result for a given observable may be used to compute that observable at a 
higher scale (i through QCD evolution. In this way the correct behavior for 
QCD radiative corrections may be incorporated using the chiral quark model 
as an initial condition for the evolution. Although radiative corrections are 
only known perturbatively, there is a chance that the matching scale, //q, 
makes perturbation theory meaningful. This a weak point which can only be 
addressed by computing higher order corrections to the model and higher 
order corrections to the QCD evolution. Nevertheless, it is encouraging 
that three different determinations of the scale based on matching to LO 
perturbative QCD evolution and available experimental data yield within 
uncertainties the scale = 320MeV for Aqcd = 225MeV. Actually, using 
this low scale the description of the distribution functions and distribution 
amplitudes agrees remarkably well with phenomenological analysis. 
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